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Resumen 


Se han evaluado, para la península de Hicacos, Matanzas, Cuba, los modelos de 
CERC (1984) y Kamphuis (2000), los cuales relacionan las tasas de transporte de 
sedimento con el flujo de energía, y el modelo de Bijker (1971), que relaciona las 
velocidades de la corriente litoral y los volúmenes de sedimento disponibles para 
el transporte. Estos han sido comparados con cubicaciones de levantamientos 
batimétricos realizados en el área. Posteriormente se comparan los resultados con 
los obtenidos mediante el software UNIBEST. 
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Introducción 


Para pronosticar el efecto de cualquier solución 
que se lleve a cabo en playas de arena es de vital 
importancia determinar la tasa de transporte de 
sedimentos. Entre los métodos más frecuentes 
para este propósito se encuentran: a) la adopción 
de tasas conocidas de áreas vecinas; b) la 
cubicación de levantamientos batimétricos de 
la playa; c) la medición indirecta por cuantifica- 
ción de volúmenes de sedimentos aportados 
por las fuentes de material y volúmenes perdi- 
dos por otros conceptos dentro de un balance 
sedimentológico de la costa; d) el cálculo directo 
mediante la estimación de las velocidades de la 
corriente litoral y la determinación de volúme- 
nes de sedimento disponible para el transporte, y 
e) el cálculo directo mediante el uso de ecuaciones 
de correlación entre transporte de sedimentos y 
el flujo de energía de los oleajes predominantes, 
mediante la aplicación de modelos numéricos de 
evolución. 

Teniendo en cuenta este último, en el 
presente trabajo se aplican las formulaciones del 


CERC (1984), Bijker (1971) y Kamphuis (2000). 
La diferencia entre ellos es la forma en que 
cada uno toma en cuenta los parámetros que 
intervienen en el proceso. Por ejemplo, con los 
modelos del CERC y Kamphuis se determina 
el transporte de sedimentos total en la zona de 
rompientes, mientras que el modelo de Bijker 
ofrece la distribución transversal del transporte 
de sedimentos longitudinal dentro y fuera de la 
zona de rompientes. 


Materiales y métodos 


Para la realización de las pruebas se definió un 
escenario de clima marítimo que consiste en las 
olas en aguas profundas generadas por el viento, 
provenientes del océano Atlántico y el Golfo de 
México. En estos cálculos se toman en cuenta 
cuatro direcciones del viento: oeste, noroeste, 
norte y noreste. 

El Instituto de Oceanología de Cuba es res- 
ponsable de investigaciones anteriores rela- 
cionadas con la erosión en la península de 
Hicacos. Plantea una tendencia erosiva en toda 
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la playa a lo largo de la península por casi 
treinta años, estimando una tasa de erosión de 
la playa mediante la comparación de perfiles a 
lo largo de los años. La conclusión a la que se 
arribó es que hay una pérdida anual de 50 000 
m' de arena por año. 


Método del CERC (1984) 


Fue desarrollado por el Centro de Inves- 
tigaciones de Ingeniería Costera (CERC, según 
sus siglas en inglés) del Cuerpo de Ingenieros 
del Ejército de Estados Unidos de América. El 
modelo está basado en la idea de que la ola que 
se aproxima desde el océano genera transporte 
de sedimentos Esto indica 
una correlación entre la tasa de transporte 


longitudinal. 


longitudinal S, y la componente longitudinal 
del flujo de energía en el borde exterior de la 
zona de rompientes. 

El modelo del CERC sólo necesita como 
entrada de datos las características de la ola 
entrante y el perfil de la costa. Los resultados de 
la fórmula del CERC dependen de coeficientes 
empíricos, tales como el coeficiente A y el índice 
de rompientes y. Este modelo es válido sólo 
para playas relativamente largas y estrechas, 
donde las diferencias a lo largo de la costa en la 
altura de la ola rompiente sean pequeñas. Más 
aún, puede ser aplicado a playas con arenas con 
diámetro medio en el rango de 0.2 a 1.0 mm. 


Fórmula usada por CERC 


Existen muchas versiones diferentes de esta 
fórmula; en este trabajo se presenta según la 
ecuación (1). 


S,= A: Hiy*co: K?-sin(0,)-cos(p,) (1 


Donde A es un coeficiente; S, el transporte 
de arena en m?/s; Cy la celeridad de la ola en 
aguas profundas en m/s; bb, el ángulo entre los 
contornos de la profundidad y la cresta de la 
ola en la zona de rompientes; H,, la altura de 
ola significativa en aguas profundas en m, y K, 
es el coeficiente de refracción. 


Suposiciones y parámetros 


El rango del coeficiente A, dependiendo del 
tipo de playa, oscila entre A =0.028-0.070 m*/s. 
El Manual de Protección de la Orilla (SPM 
según sus siglas en inglés) recomienda un valor 
de 0.05, mientras que el valor más usado en 
Holanda es de 0.040. En los cálculos se usará 
un valor de 0.030. 


Método de Kamphuis (2000) 


El modelo de Kamphuis es una variante de 
la fórmula del CERC obtenida en un estudio 
canadiense de sedimentos costeros. Kamphuis 
(1991) llegó a una expresión que incluye los 
efectos del periodo de la ola (o esbeltez de la 
ola), la pendiente de la playa y el diámetro del 
grano. 

La fórmula de Kamphuis es válida para 
playas planas, sin complicaciones en el sistema 
de barreras de rompientes, en donde no 
existan corrientes de marea, con contornos de 
profundidad rectos, sin presencia de campos 
de groins ni rompeolas fuera de la zona hacia 
fuera de la orilla. 


Fórmulas usadas para Kamphuis 


Un análisis dimensional de la tasa de sedimen- 
to Q. en la turbulenta zona de rompientes 
conduce a la siguiente expresión general: 


Q 


E -1.25 E 0.25 

$ o 3 b 0.75 b 

o (es) Mm, (ls 
Ash 

dd o) 

¡me*(o -0Ly) 


En la cual Q es el transporte de sedimentos 
longitudinal a la costa en kg/s; p, la densidad 
del agua del mar en kg/m; H,, la altura de la 
ola en la zona de rompiente en m; L, el periodo 
pico del espectro en s; L,, la longitud de la ola 
en aguas profundas en m; m, la pendiente de 


la playa; D,,, el diámetro medio del grano en 
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m, y Ol, es el ángulo de la ola en el punto de 
rompiente en grados. Para obtener el transporte 
de sedimento en unidades métricas (m*/s) se 
debe dividir Q por la densidad del sedimento y 
la fracción de volumen ocupada por el material 
(I-porosidad): 


Q 
de == 
(1—p)p, 


Donde S es el transporte de sedimento 
en m?/s; p, la porosidad en tanto por 1, y p, 
es la densidad del sedimento en kg/ m?. Esto 
conduce a la siguiente fórmula: 


¿13x10% pH (m0 Y" 
(—p)p, E Lo 


(3) 


0.25 
7073 (255) sen0S(24b) 


Suposiciones y parámetros 


Según Kamphuis et al. (1986), el valor del 
coeficiente de Kamphuis fue determinado en 
1.3-10* o 41 000 para resultados por año. De 
acuerdo con Kamphuis (1991), este valor debe 
ser 64 000. Schoonees (1996) analizó más datos 
y propone un valor de 1.6*10* (o 50 000) para 
condiciones de oleaje suave, y/o material 
grueso. Éste debe ser el valor correcto para 
calcular el transporte de sedimentos potencial 
por año. Para sitios expuestos, el mejor valor 
es 65 000. 


Método de Bijker (1971) 


Este modelo calcula el transporte de sedi- 
mentos como una función del campo de 
oleaje dado y de la corriente longitudinal 
dada, independiente del modelo original (olas 
inducidas por corriente O por mareas). La 
fórmula es capaz de predecir la distribución 
del transporte de sedimentos sobre el perfil, así 
como el transporte total. 


La componente de carga de fondo fue 
adaptada de una fórmula existente para la 
carga de fondo bajo circunstancias de río. 
Bijker dividió la expresión en un parámetro de 
remoción y otro de transporte, e introduce la 
influencia del oleaje para modificar la tensión 
de cortante del fondo, T, y convertirla en el 
parámetro de remoción T,.. 

Este método tiene como ventajas que se 
resuelve cualquier problema que también se 
pueda con la fórmula de CERC; es adaptable a 
cualquier condición de corriente; el concepto de 
la fórmula de Bijker —el ajuste de la tensión de 
cortante tenida en cuenta para las olas— puede 
tener una aplicación mucho más universal; 
alternativamente, la velocidad, V, incluida 
en la fórmula de Bijker, puede ser conducida 
por cualquier combinación de fuerzas y sujeta 
a toda clase de influencia local, y la fórmula 
de Bijker puede ser usada para predecir la 
sedimentación en un canal estrecho, en el cual 
no hay olas rompientes (Van der Velden, 2000). 


Fórmulas usadas para Bijker 


El modelo de Bijker está basado en el concepto 
de que el transporte de sedimentos “5” es igual 
al producto de la velocidad del flujo *V” y la 
concentración de sedimentos “C”: 
S=5),+S, (4) 
En la cual S, es el transporte de arrastre 
del fondo y S, es el transporte de partículas 
suspendidas. 


La formulación del transporte de fondo S, 
es: 


, ED VE o -0.27:A-Dsy"P-8 | (5) 
á o UT 


En la cual S, es el transporte de fondo en 
m?/sm; b, un coeficiente; D, el diámetro de la 
partícula en m; V, la velocidad local en m/s; r, 
la rugosidad de fondo en m, g, la aceleración de 
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la gravedad en m/s?; y, el factor de ondulación 
= (C/C90); C, el coeficiente de Chézy, y A es 
la densidad relativa aparente del material 
de fondo. La tensión de radiación de fondo 
promediada en el tiempo debido a oleaje y 
corriente “7, ' es: 


113 
07 loz (5-2) | (6) 


Tew = 


ES 


En donde £¿ es el parámetro de Bijker, 
=C PE ; fw es el factor de fricción de Johnson, 
8 


= expl-5.977 + 5.213(a0/r)-0.194]; C es el 
coeficiente de Chézy en Vm/s, = 18 log(12 
h/r); C90 es el coeficiente de Chézy para el 
diámetro 90 en Vm/s, = 18 log(12 h/D90); h es 


la profundidad local en m, uo es la velocidad 


orbital en m/s, =a0 21 ; no es el desplazamiento 


máximo en el fondo en m, = ETICO! 


Para determinar el transporte suspendido 
conociendo el transporte de fondo es necesario 
determinar el factor de las integrales de Eins- 
tein, Q, como se aprecia en la ecuación (8), para 
lo cual se utilizan los valores de las integrales 
de Einstein, donde se debe calcular el número 
de Rouse, z* y la relación r/h. 


Z'= (7) 


Donde x es la constante de von Karman = 0.4 
y V*es la velocidad de esfuerzo cortante 
[m/s]. 

Ésta es adaptada, teniendo en cuenta la 
influencia de las olas. =1* ¿q = es 


El sedimento suspendido está directamente 
relacionado con la carga de transporte de fondo 
de la siguiente manera: 


S, =1.83-Q-S, (8) 


El transporte de sedimentos en toda la 
sección se determina mediante: 


S =S»(1+1.83-Q- S,) (9) 


Distribución de velocidad en la zona de 
rompiente 


Si se analiza una columna de agua existirán 
fuerzas que provocan el movimiento de las 
aguas y fuerzas que se oponen a éste. El 
equilibrio de estas fuerzas en la dirección 
a lo largo de la costa provoca una corriente 
constante en esa dirección. 

Como fuerzas que provocan el movimiento 
se encuentran las componentes de la ten- 
sión de cortante, las fuerzas del viento y 
las provocadas por la marea. Como fuerzas 
resistentes turbulencia 
(mecanismo de difusión horizontal) y la fuerza 
de fricción del fondo. 

Una predicción de la distribución de 
velocidad usada en el modelo original de Bijker 
(desde ahora se referirá como Bijker, original) 


se encuentran la 


está basada en el equilibrio entre la fricción del 
fondo y tensión de cortante. Estableciendo este 
equilibrio se podrán predecir los resultados de 
la velocidad en un perfil lineal perpendicular 
a la costa. 


v= 5 sen(po) C 
8/2 co  ffw 


v/g-h-m (10) 


Si se acepta este planteamiento, entonces 
la velocidad de la corriente longitudinal es 
una función lineal de la profundidad del agua 
dentro de la zona de rompiente, ocurriendo el 
máximo de velocidad en el borde externo de la 
zona de rotura. Fuera de la zona de rotura la 
velocidad es cero. 

Longuet-Higgins (1971) ha predicho teóri- 
camente la distribución de velocidades, 


incluyendo la fuerza turbulenta en un 
equilibrio dinámico junto con la tensión 
de cortante y la fricción. Según la teoría de 
Longuet-Higgins (1971), los efectos de fricción 
lateral de la corriente longitudinal producen 


una distribución transversal no uniforme a la 
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velocidad del flujo, con un máximo entorno a 
la línea de rompiente (Medina, 1997). 


+  Normalizar la abcisa, Y = y/y, donde 
y, indica, con la notación habitual, la 
coordenada del punto de rotura (en 
este caso, la distancia a la línea de la 
orilla). 

e Calcular el parámetro P, definido como 


Tm:N, donde N es una constante 
20:-C 


adimensional que ha sido acotada 
superiormente por N =0.010 por Galvin 
y Eagleson (1965), por lo que puede 
tomarse este valor a falta de uno mejor; 
Q es una constante empírica, cuyo 
valor se sitúa en torno a 0.4; C es un 
coeficiente representativo de la fricción 
del fondo, su valor se sitúa en torno a 
0.010. 
e Calcular los parámetros: 


A NAO 1 
=|-— + + , = , 
4) V16 1-2.5P 
ollas MEL y 
pl- p2 p1- p2 


e Determinar la distribución transversal 
adimensional de la velocidad: 


V* (y) =A-Y+B1-Y” a an) 
V"(4)=8B2+ Y l<Y<o 


donde V*(y) esla velocidad litoral adimensional 

en el punto y [-] = V/Vb. 

e Determinar la distribución transversal de 
la velocidad: 


V(y) = V*(y) - Vb (12) 


Donde Vb es la velocidad litoral en la línea 
de rompiente en m/s. 

Si en realidad el oleaje no es regular, la 
altura de la ola puede variar y el borde exterior 


de la zona de rompientes cambiar de posición. 
Este efecto trae como consecuencia que el perfil 
de velocidades se haga más ancho y con menor 
pico en su forma. 

Estas dos aproximaciones de distribución 
de la velocidad de la corriente longitudinal 
se emplean en los cálculos; la primera, la 
expresión simple, es usada en Bijker (original) 
y la segunda, con la distribución de Longuet- 
Higgins, en Bijker (LH). 


UNIBEST 


UNIBEST es un softivare creado por Deltares, 
que ha sido empleado en la modelación de 
varios proyectos costeros en Holanda y otras 
partes del mundo con muy buenos resultados. 

Para calcular el transporte de sedimentos 
paralelo a la costa puede utilizarse el módulo 
UNIBEST LT (Long shore transport). Puede 
calcularla cantidad de transporte de sedimentos 
debido a condiciones de corriente y oleaje, y 
los datos que provienen de estos cálculos son 
usados en UNIBEST CL (Coast line development) 
para predecir cómo se desarrollará la línea de 
costa en un periodo de tiempo dado y cómo 
ello afecta el transporte local. 


Caso de estudio 


La península de Hicacos, situada en la pro- 
vincia de Matanzas, Cuba, presenta serios 
problemas de erosión. Tal es así, que la 
playa de Varadero, que abarca los 20 km 
de la costa norte de la misma, está siendo 
constantemente monitoreada y atendida para 
no perder tan importante recurso natural y 
turístico, el más importante del país. Esta 
playa presenta un grupo de características 
que la hacen representativa no sólo del resto 
de las playas de Cuba, sino también de la 
región del Caribe, algunas de ellas son la 
dirección del clima marítimo, destacándose 
la dirección O, NO, N y NE; la composición 
del sedimento, predominando los restos 
de organismos vivos; la similitud entre 
variables climatológicas, como temperatura, 
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radiación solar y precipitaciones, entre otras, 
correspondientes a un clima tropical; y la 
importancia socioeconómica de la misma, 
utilizada como centro turístico y zona de 
asentamiento poblacional. 


Oleaje de entrada 


De acuerdo con la literatura, los cálculos con 
estos métodos se pueden hacer por dos vías 
diferentes de entrada de datos del oleaje: la 
primera tomando una ola que responda al 
promedio ponderado para cada dirección, O, 
NO, N y NE, y la segunda, para todo el clima 
de oleaje de estas direcciones. 


e Primero, una ola media ponderada por 
dirección del viento es tomada en cuenta 
para facilitar los cálculos. Para esto, los 
datos de la dirección, altura y periodo de 
la ola y su probabilidad de ocurrencia se 
tomaron del Global Wave Statistics (GWS), 
zona 32. La altura de ola significativa 
media ponderada se calcula multiplicando 
la altura de cada ola por su probabilidad 
de ocurrencia. Esto es hecho para todas las 
alturas de ola y la suma es el valor para 
esa dirección. Lo mismo es hecho para 
determinar los periodos de la ola. Los 
resultados se muestran en el cuadro 1. En 
los cálculos, cada ola causa un transporte 
proporcional a la ocurrencia de la dirección 
representativa. 

e Segundo, se toma en cuenta todo el 
clima de oleaje que afecta la península de 
Hicacos, según el Global Wave Statistics, 
zona 32, en la cual se puede obtener cada 


altura de ola posible y su probabilidad 
de ocurrencia. Para cada rango de alturas 
de ola se ha calculado un periodo medio 
multiplicando el periodo de cada ola por 
su por ciento de ocurrencia. Esto significa 
un clima de oleaje de 37 olas significativas 
en aguas profundas. Cada una de estas olas 
específicas produce una cierta cantidad de 
transporte de sedimentos. 


Sólo hay una corriente de marea débil 
presente en la península. Si la velocidad de 
esta corriente es despreciada, esto no tendrá un 
gran impacto en el sistema. 

El valor del índice de rompiente que se toma 
es de 0.8. 


Características de la arena 


A partir de los estudios realizados por el 
Instituto de Oceanología (Juanes, 1997), se 
determina el resto de los parámetros usados en 
estos cálculos relacionados con el sedimento y 
que se muestran en el cuadro 2. 

La fuente principal de los sedimentos son 
las algas Halimeda algae. 


Batimetría y transporte de sedimentos 


En Varadero existe una batería de perfiles 
perpendiculares a la playa que se muestran 
en la figura 1. Estos perfiles cubren 
aproximadamente hasta 50 m al interior de la 
playa, tal y como se observa en las figuras 2a y 
2b. Según Juanes (1997), un estudio realizado 
a tres perfiles de la playa sirvió para hacer una 


primera aproximación del ritmo de pérdida 


Cuadro 1. Olas promedio obtenidas del GWS. 


Dirección Ángulo de aproximación 0, () | Altura de ola significativa media Hs (m) | Periodo medio, T (s) 
Noreste -75 1.51 5.10 
Norte -30 1.71 5:92 
Noroeste 15 1.74 5.19 
Oeste 60 SS 4.60 
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Cuadro 2. Parámetros usados en los cálculos del modelo de Bijker. 


Símbolo Parámetro Valor Unidad 
Pesa Densidad de la arena 2700 (kg/m?) 
pa Densidad del agua 1025 (kg/m?) 

De, Diámetro del grano medio 0.26 (mm) 
¡Di Diámetro 90% del grano 0.40 (mm) 
r Coeficiente de rugosidad 0.050 (m) 
w, Velocidad de caída del grano 0.04550 (m/s) 
IV TH 
TI 
I 
xv XIV Perfiles : Puntos de muestreo 
topográficos de areña 
L Punta Molas 
> TL W de Punta Molas 
ys TIL Guardafronteras 
eo IV. Rincón Francés 
Esc. 1:150 000 e Me “elas 
Oo vVL Las Américas 
Oo VII Villa Cuba 
VIL  Caney 
IX. Villa Los Cocos 
a) X. Calle 16 
XL Villa Tortuga 
XIL Kawuama 
XUL Punta Blanca 
To XIV. Hotel Oasis 


XV. La Conchita 


Figura 1. Distribución de los perfiles y puntos de muestreo a lo largo de la península de Hicacos (Juanes, 1997). 


general, considerándose una longitud de playa 
de 20 km, obteniéndose que la pérdida de 
volumen está en el orden de los 40 000 a 50 000 
m'/año. 


Perfil promedio 


Para simplificar los cálculos se ha seleccionado 
un perfil medio para toda la península. Este 


perfil es tomado cerca del Hotel Bella Costa, el 
cual se muestra en la figura 3. Se ha agregado 
una línea recta en este perfil para simplificar 
los cálculos. Con esto se asume que esta línea 
es representativa para todo perfil de la playa 
en toda la península. La línea de tendencia 
tiene una pendiente de 1:60. Es decir, por 
cada metro de profundidad se ha medido una 
distancia de 60 metros. 
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a) Perfiles para la sección PN-18 
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b) Perfiles para la sección PN-24 


Figura 2. Perfiles medidos en la península. Fuente: Instituto de Oceanología, Cuba. 


gua, vol. Il, núm. 3, julio-septiembre de 2011 


Ciencias del A 


d vJB0]0u92], 


Córdova-López y Torres-Hugues, Modelo matemático para la determinación del transporte longitudinal para playas del Caribe 


0 200 400 


Profundidad (m) 


15 


18 


Distancia (m) 


600 800 1000 1200 


Figura 3. Perfil medio usado para los cálculos 


Resultados y discusión 
Ancho de la zona de rompientes 


En primer lugar, se determina el ancho de la 
zona de rompientes y valores de algunas de 
las características más importantes en el resto 
de los cálculos para las olas representativas de 


cada dirección, tal como aparece en el cuadro 
3. 


Resultados de cálculo con el modelo del CERC 


Teniendo en cuenta el rango de aplicación 
del modelo y las características de la playa en 
estudio, se puede afirmar que este modelo es 
aplicable a esta playa. 

Usando el modelo del CERC se han deter- 
minado las siguientes cantidades de transporte 
de sedimentos con los datos mostrados en el 
cuadro 3 para las cuatro olas, una por cada 


dirección, y con todo el clima de oleaje. El 
resultado del cálculo se tiene en el cuadro 4. 


Resultados de cálculo con el modelo de 
Kamphuis 


En el área de interés (península de Hicacos) 
existe alguna influencia de la marea, pero no 
es mucha realmente (la diferencia de mareas 
es de 0.20 m). Hay rompeolas en la salida 
del Canal de Paso Malo y los contornos de 
profundidad difieren a lo largo de la penínsu- 
la; pero tomando esto en consideración, el 
modelo de Kamphuis aún puede ser usado 
para una estimación de la cantidad potencial 
de transporte de sedimentos. Para estos 
cálculos, el coeficiente de Kamphuis ha sido 
determinado como 1.3-10* m/s o 41 000 
m/a, porque en la playa de Varadero, las 
condiciones de oleaje no son suaves y no hay 
presente material de fondo grueso. 


Cuadro 3. Características fundamentales en la línea de rompiente para las cuatro olas. 


Dirección 0, 4) Hs (m) T (s) ZR (m) Hb (m) hb (m) Vb (m/s) 
Noreste -75 1.51 5.10 26 0.91 1.14 0.87 
Norte -30 11ZAL 5:32 103 1.67 2.09 1.86 
Noroeste 15 1.74 5.19 106 1.75 2.18 2.00 
Oeste 60 159 4.60 Ya 1.05 1,51 Ts 
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Usando este modelo se calcularon las si- 
guientes cantidades de transporte de sedi- 
mentos con los datos del cuadro 3 para cuatro 
olas, una para cada dirección, y para todo el 
clima de oleaje. Los resultados se presentan en 
el cuadro 4. 


Resultados de cálculo con el modelo de Bijker 


Para los cálculos por este modelo se toman los 
valores recomendados para el coeficiente b = 1 
para la zona fuera de la zona de rompiente, y 5 
para el interior de la zona de rompiente. 

Los cálculos se hicieron tanto con la 
distribución de la velocidad con la teoría simple 
como con la distribución de la velocidad de 
Longuet-Higgins. 

Usando el modelo de Bijker se calculan 
las siguientes cantidades de transporte de 
sedimentos con los datos del cuadro 3 para las 
cuatro olas, una para cada dirección, y para 
todo el clima de oleaje. Los resultados se tienen 
en el cuadro 4. 


A modo de resumen se presentan los 
resultados de los cálculos en el cuadro 5. 
Observando estos, se puede decir que: 


e El signo negativo indica el sentido del 

transporte, tomándose como convenio 
que el transporte hacia el oeste es positivo 
y hacia el este es negativo. En este caso 
se puede apreciar que el transporte neto 
total es positivo, por lo que el transporte 
neto es en la dirección oeste, es decir, 
desde Cayo Francés hacia el Canal de Paso 
Malo, comportándose según los estudios 


realizados por el Instituto de Oceanología. 


Si se observan los valores del transporte 
de sedimento se pueden notar cuestiones 
interesantes: 


e La diferencia entre los valores calculados 
con las cuatro olas representativas y todo 
el clima de oleaje se debe a la dependencia 
no lineal de la altura de ola rompiente en 


Cuadro 4. Resultados del cálculo de la tasa de transporte. 


Dirección P (%) a Ss Vara A cmphuis Ss ora A nieto Ss Teixertor s A ikerttÉn Ss A etom 
del oleaje (m*/año) (m*/año) (m*/año) (m'/año) (m'/año) (m*/año) (m'/año) (m*/año) 
Noreste 18.19 312 000 590 000 46 500 66 670 112 000 218 100 55 000 107 100 
Norte 9.82 652 000 1001 000 24 000 49 380 179 000 293 800 61 000 100 100 
Noroeste 4.93 -154 000 -301 000 -13 900 -24 500 -50 000 -82 700 -28 000 46 300 
Oeste 3:29) -86 000 -186 000 -5 100 -10 240 -37 000 -78 000 -7 000 -14 800 
Total 724 000 1 104 000 51 500 81300 204 000 351 200 81 000 146 100 
Donde P : por ciento de ocurrencia; 54... : tasa total de transporte de sedimentos usando cuatro olas representativas mediante el método de CERC; 
ST espe : tasa total de transporte de sedimentos usando todo el clima marítimo mediante el método de CERC; 54, .,,,, : tasa total de transporte 


de sedimentos usando cuatro olas representativas mediante el método de Kamphuis; ST, 


todo el clima marítimo mediante el método de Kamphuis; S4 
mediante el método de Bijker(OR); ST 


Bijker(OR) * 


Bijker(OR); S4 : tasa total de transporte de sedimentos usando cuatro olas representativas mediante el método de Bijker(LH); ST 


Bijker(LH) * 


Bijker(OR) * 
tasa total de transporte de sedimentos usando todo el clima marítimo mediante el método de 


: tasa total de transporte de sedimentos usando 


Kamphuis 


tasa total de transporte de sedimentos usando cuatro olas representativas 


Bijker (LH) * 


tasa total de transporte de sedimentos usando todo el clima marítimo mediante el método de Bijker(LH). 


Cuadro 5. Resumen de los resultados de los cálculos. 


Capacidad de transporte de Bijker Bijker 
sedimentos (m//a) CERC (original) (Longuet-Higgins) Kamphuis 
Cuatro olas representativas 724 000 204 000 81 000 51 500 
Todo el clima de oleaje 1 104 000 351 200 146 100 81 300 
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la ecuación de transporte de sedimentos. 
Esto conduce a resultados mayores. Si, por 
ejemplo, las olas siguen una distribución 
tipo Rayleigh y las olas tienen periodo 
uniforme, la tasa de transporte calculada 
usando la distribución de alturas de olas 
más 
grande que los calculados usando sólo una 


es aproximadamente 1.53 veces 


altura de ola media ponderada. Esto explica 
los resultados más grandes. Cuando se 
usa todo el clima de oleaje, las olas están 
representadas con más exactitud, causando 
mayores capacidades de transporte de 
sedimentos y, por tanto, mayores tasas de 
erosión. 

e En cuanto a las tasas de transporte, se 
observó que los cálculos realizados con la 
formulación del CERC daban valores muy 
elevados, mientras que la formulación de 
Kamphuis daba valores muy pequeños. 
Estos dos métodos son dos métodos 
integrales que no describen el transporte 
de sedimentos a partir de la descripción 
detallada de los fenómenos que tienen 
lugar. Son fórmulas generales calibradas 
con mediciones de laboratorio y prototipo, 
por lo que la universalización de las mismas 
puede brindar resultados no satisfactorios, 
tal y como ocurre en el caso de la playa de 
Varadero. 


El modelo de Bijker usa la tensión de corte 
del fondo para calcular la distribución en el 
perfil, mientras que los modelos de CERC y 
Kamphuis emplean la altura de la ola entrante 


y su dirección, usando un modelo de energía 
de más, mostrando que no hay diferencia en 
la cantidad de transporte de sedimentos a 
través de una unidad de ancho en una sección 
transversal. 

Los resultados del transporte desedimentos 
longitudinal neto para los primeros cincuenta 
metros de playa se observan en el cuadro 
6. Es de destacar que las estimaciones del 
transporte de sedimentos realizadas por el 
Instituto de Oceanología, a partir de los per- 
files batimétricos existentes en Varadero, 
plantean un desplazamiento de la arena entre 
los 40 000 y los 50 000 m*/año, el cual puede 
comprender otros fenómenos, tales como el 
transporte transversal, producto de eventos 
extremos, que pueden conducir parte de la 
arena fuera del sistema y por tanto no retorna 
al origen. Observando las figuras 2a y 2b se 
puede apreciar que aproximadamente la mitad 
de los perfiles corresponde a la parte seca de 
la playa, por lo que se puede estimar que el 
transporte de sedimentos longitudinal es un 
70% del transporte total, por tanto debe estar 
en el orden de los 28 000 a 35 000 m*/año. En 
el cuadro 6 se muestran los resultados para el 
resto de las formulaciones, pero integradas 
solamente para los primeros cincuenta metros 
de playa, lo cual fue realizado por Trujillo 
(2004). De estos resultados se observa que el 
método de Bijker, utilizando la distribución de 
corrientes longitudinales original (propuesta 
por Bijker) y la formulación propuesta por 
Longuet-Higgins (1971) para el cálculo de las 
corrientes longitudinales generadas por el 


Cuadro 6. Transporte de sedimentos para los primeros cincuenta metros (m*/año). 


Dirección CERC Kamphuis Bijker (OR) Bijker (LH) 
NE 183 600 14 100 35 300 27 000 
N 117 000 8 300 7400 11 200 
NO -28 400 -2 600 -3 200 -5 600 


O -40 500 -2 400 -6 300 -3 400 
Total 231 700 17 400 33 200 29 200 
Oceanología entre 28 000 y 35 000 


Tecnol 


Ciencias del Agua, vol. IL, núm. 3, julio-septiembre de 2011 


0 gía y 


gua, vol. Il, núm. 3, julio-septiembre de 2011 


Ciencias del A 


d vJB0]0u92], 


Córdova-López y Torres-Hugues, Modelo matemático para la determinación del transporte longitudinal para playas del Caribe 


oleaje, brinda un resultado muy satisfactorio 
(33 200 y 29 200 m*/año, respectivamente), en 
comparación con las estimaciones del Instituto 
de Oceanología; mientras que la formulación 
del CERC brinda un resultado muy grande (231 
700 m*/año) y la formulación de Kamphuis da 
un resultado pequeño (17 400 m*/año). 


Comparación de un modelo simplificado en 
UNIBEST con los resultados de los cálculos 
anteriores 


A continuación se lleva a cabo una comparación 
entre los resultados con el mismo modelo 
simplificado arrojados por UNIBEST LT y 
los resultados estimados anteriormente, para 
comprobar si el modelo UNIBEST arroja 
resultados similares a los cálculos con Bijker 
(LH). 

Como se había visto, en el primer caso se 
ha tomado sólo una ola representativa de cada 
dirección; en el segundo caso se ha corrido 
todo el clima de oleaje. Se ha usado el mismo 
perfil. El borde dinámico se ha puesto a una 
profundidad de 11 m (672 m fuera de la costa). 

Con la aproximación de Bijker se calculó 
la cantidad y dirección del transporte de 
sedimentos. Para ello se utilizaron los mismos 
parámetros. Los resultados de la corrida con 
LT en UNIBEST se presentan en el cuadro 7. 
El modelo UNIBEST arroja un valor negativo 
en la dirección oeste, mientras que el estimado 
dio un resultado positivo para esta dirección. 
Por tanto, el signo de este estimado ha sido 
cambiado en el cuadro 7. 


Observando los resultados del cuadro 
7 relacionados con la determinación de la 
capacidad de transporte de sedimentos, se 
puede notar que los resultados están en el 
mismo orden de magnitud entre los valores 
obtenidos mediante el modelo de Bijker (LH) y 
el modelo UNIBEST tanto para las cuatro olas 
representativas de las direcciones como para 
todo el clima de oleaje marítimo. 

Esto, sin duda, guarda estrecha relación con 
la configuración y los valores de la distribución 
de la velocidad de la corriente longitudinal 
mostrada en la figura 4. De la observación de 
esta figura se puede señalar que ambos modelos 
presentan aproximadamente la misma área bajo 
la curva, y que la velocidad según Longuet- 
Higgins presenta valores mayores, de ahí que 
también se produzcan mayores valores en la 
tasa de transporte de sedimentos longitudinal 
que los obtenidos con el UNIBEST, pero en el 
mismo rango de magnitud. 

Si se compara la distribución de velocidad 
encontrada en UNIBEST con la de los cálculos, 
se pueden apreciar diferencias entre las usadas 
en Bijker (original), Bijker (Longuet-Higgins) 
y Bijker (UNIBEST). Dicha comparación se 
muestra en la figura 4. 

El modelo UNIBEST usa para describir la 
distribución de la corriente longitudinal en la 
aproximación de Bijker la ecuación básica de 
momentum (ecuación (13)): 


d dh 8 
—= 8, +pgd— + VIV, =0 13 
de Y Pg dy p C? | h A (13) 


Cuadro 7. Comparación del UNIBEST con los resultados anteriores. 


A os con NEBESn 
Bijker (LH) (m*/y) Bijker (m//y) 

NE -55 000 -37 600 

N -61 000 -68 500 

NO 28 000 +22 200 

(0) 7000 +10 100 

Total -81 000 -84 000 

Todo el clima marítimo -146 100 -121 400 
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Distribución de la velocidad de la corriente (NE) 


1.800 


1.200 


0.600 


Velocidad (m/s) 


0.000 


—— V (orig.) 
—S- V (L-H) 


—+- V (UNIBEST) 


Distancia desde la costa (m) 


Figura 4. Distribución de la velocidad de la corriente en la dirección NE. 


En la cual, la segunda parte (dh, / dy) invo- 
lucra a la marea. Pero como se puede observar 
de la figura4, el área bajo el gráfico para Longuet- 
Higgins y UNIBEST son aproximadamente 
las mismas. De esto se puede concluir que se 
toma en cuenta la misma cantidad potencial de 
transporte de sedimentos. 


Conclusiones 


El modelo Bijker (LH) es capaz de describir 
con buen grado de aceptación el mecanismo 
del transporte de sedimentos longitudinal, 
por lo que se recomienda que este modelo sea 
utilizado en el resto de las playas de la región, 
ajustándolo a las características particulares de 
cada playa. 


Recibido: 21/10/09 
Aprobado: 30/09/10 
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Abstract 


CÓRDOVA-LÓPEZ, L.E € TORRES-HUGUES, R. Mathematical model to determine 
alongshore transport rate to Caribbean beaches. Water Technology and Sciences, formerly 
Hydraulic engineering in Mexico (in Spanish). Vol. IL, No. 3, July-September, 2011, pp. 
127-140. 


This article evaluates CERC (1984) and Kamphuis (2000) models for the Hicacos Peninsula, 
Matanzas, Cuba, which relate sedimentation transport rates with energy flow. The Bijker 
(1971) model is also evaluated, relating alongshore current velocity with sediment volumes 
available to be transported. These are compared with volumetric/bathymetric measurements 
conducted in the area. The results are then compared with those obtained using UNIBEST 
software. 


Keywords: transport, longitudinal, sediment, Caribbean, beach, model, current. 
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